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Abstract 

We introduce an xy generalization of the frustrated Ising model on a tri- 
angular lattice. The presence of continuous degrees of freedom stabilizes a 
finite-temperature spin state with power-law discrete spin correlations and an 
extensive zero-point entropy. In this phase, the unquenched degrees of free- 
dom can be described by a fluctuating surface with logarithmic height correla- 
tions. Finite-size Monte Carlo simulations have been used to characterize the 
exponents of the transition and the dynamics of the low-temperature phase. 
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I. INTRODUCTION 



The possibihty of spin systems with a residual zero-point entropy has existed since Wan- 
nier's study of the two-dimensional Ising triangular antiferromagnet;§ here the frustration 
leads to a zero-temperature (T = 0) critical point characterized by an extensive entropy 
5*0 = 0.338. The associated unquenched degrees of freedom can be described by a fluctuat- 
ing surface with logarithmic height correlations;!'! in discrete spin models the weak nature 
of the long-range entropic interactions between height defects forbids the existence of this 
critical phase at T 7^ 0. Can a phase transition into a state with a zero-point entropy exist 
at finite temperatures? This issue was first raised by AndersonS in his study of spinels, 
where he noted that a [111] projection of this structure corresponds to a planar kagome net. 
Anderson suggested that the ground-state manifold of Heisenberg spins on such frustrated 
lattices would be highly degenerate, and thus that long-range magnetic order would be sup- 
pressed due to the residual entropy at finite temperatures.! The resulting phase would have 
power-law discrete spin correlations at finite temperatures in contrast to the exponential and 
constant Ising correlations associated with a paramagnet and an ordered phase; it is thus 
natural to call this state a "classical spin liquid" . 

There now exists an experimental realization of Anderson's model: the magnetoplumbite 
SrCrs-xGai+^OiQ {SCGO{x)), whose magnetic properties are attributed to planes of 
antiferromagnetically-coupled chromium atoms on a kagome lattice.! Though inelastic neu- 
tron scattering! and specific heati measurements on SCGO{x) are consistent with antifer- 
romagnetic order, the observed spin correlation length is only a few lattice spacings.! There 
exist many other anomalous magnetic materials that, like SCGO{x), do not display conven- 
tional long-range spin order.!'! Most of these experimental systems have physical structures 
that impose strong constraints on the low-temperature configurations of their continuous 
spins, allowing them to be identified with the ground-states of an associated discrete model 
on the same lattice. Following the spirit of Anderson's original suggestion,! we show that an 
interplay between continuous and discrete degrees of freedom can enhance the defect-defect 
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interaction in the associated fluctuating surface, tlius stabilizing a finite-temperature critical 
point associated with a zero-point entropy. 

II. THE DISCRETE MODEL 

The two-dimensional Ising antiferromagnet with nearest-neighbor interactions on the 
triangular lattice is the simplest frustrated spin system;0 unlike its ferromagnetic counterpart 
it does not order even at zero temperature but has power-law spin correlations in contrast to 
a conventional paramagnet.i The geometry of this lattice does not permit the minimization 
of each bond energetically and thus even at T = there exist parallel and antiparallel 
nearest-neighbor spin pairs which we denote as antiferromagnetic (AFM) and ferromagnetic 
(FM) bonds for convenience. The energy of each triangle is minimized if it contains one and 
only one ferromagnetic bond. The number of spin configurations that satisfy this plaquette 
constraint scales exponentially with the site number; the resulting T = manifold has an 
extensive entropy and can be fully explored from a given starting state by flipping spins 
with equal numbers of AFM and FM bonds (Figure 1). If a spin with four AFM bonds is 
reversed there result two "bad" triangles that violate the constraint (Figure 2). These two 
"defects" can be separated to arbitrary distance at no additional cost, despite the fact that 
in doing so all the spins along the line linking them are reversed (Figure 3). The situation 
is similar to that in the Id Ising ferromagnet where a flip of a single spin creates two "bad" 
bonds. These defects can move freely along the chain by reversing all the spins between 
them; since no more "bad" bonds appear in this process their energy cost is independent 
of their relative positions. In the triangular case the analogous separation of two "bad" 
triangles does not lead to the formation of more "bad" plaquettes, so that their energy 
cost is also finite in the thermodynamic limit. At T = the creation of these defects is 
forbidden, there is no length-scale in the problem and the spin correlations are power-law. 
However at nonzero temperatures defect formation is favored entropically, and the spin 
correlations become exponential with a correlation length determined by the defect density. 
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We contrast this situation with that of the 2d Ising ferromagnet where the minimum energy 



is suppressed until a finite temperature. 

We can characterize each spin state by its spatial configuration of antiferromagnetic 
bonds. At T = the plaquette constraint ensures that each triangle has exactly two AFM 
bonds; each spin state can thus be mapped onto a unique tiling of the plane by rhombi 
which can also be viewed as a planar projection of a surface in three-dimensions (Figure 4). 
Quantitatively this mapping! between the spins {si} and the heights {hi} associated with 
the surface is defined by the expression 



for sites i and j as shown in Figure 4. If the underlying spin configuration satisfies the 
plaquette constraint, the sum of the height differences as defined in (|l|) around any closed 
path will be exactly zero (Figure 5a) and the height is a uniquely defined function. By 
contrast, a similar sum determined along a loop enclosing a "bad" triangle (see Figure 5b ) 
results in a finite height difference which, in analogy with conventional elasticity theory, we 
denote as the Burger's vector; the underlying spin configuration determines its numerical 
value to be a multiple of six. Defects in the surface appear because the height is now a multi- 
valued function whose branching points correspond to triangles that violate the constraint. 

The surface representation provides a very convenient description of the full ground-state 
manifold; in this language the spin fiips associated with a zero energy cost correspond to 
the addition/subtraction of elementary cubes. Thus the allowed T = spin configurations 
are mapped onto a free fiuctuating surface whose behavior at finite temperatures and long 
length-scales is characterized by the free energy density 



cost associated with two defects at a separation L scales like E ^ L^-^ here defect formation 




(1) 




(2) 



where rjT is the stiffness associated with the surface tension of the membrane. Equation (^ 



leads to logarithmic height correlations of the form 
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m-h,)') = —\n'^ (3) 

TjTT a 

where a is the cutoff imposed by the lattice. Since the spins and the heights on any sites i 
and j are related by the expression 



Si ■ s 



the spin-spin correlation function at large distances is power-law 

TV 

7J-2 ( a\'^^ 

{si ■ Sj) = exp -Y((^i - hjf) = ( — 1 (5) 

which follows from the free energy density (0). For the T Ising model with nearest- 
neighbor interactions that we have considered above, the correlation function (^) is known 
from the exact solutioni'i and results in a stiffness coefficient 

Vising = ^ (6) 

for the associated free surface. In principle, a generalized spin model on the triangular lattice 
could lead to a different numerical value of rj, which is determined by the specific form of 
the spin-spin interaction. 

The free energy density (0) leads to a defect energy that scales logarithmically with 
system size at any finite temperature. The resulting free energy associated with an isolated 
defect is 

Fdefect = ^0 - 2T I (^^^ - 1 1 In L (7) 

where Eq is the core energy and the last term in (|^ is the entropy associated with the 
random location of the defect on the surface. This defect free energy is logarithmically large 
if the coefficient of the temperature in (0) is positive; for the case here where the minimum 
Burger's vector bmin = 6, the critical value for r] is 

Stt 2n 
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Binding of the surface defects^ occurs for r] > r]c^, and we note that rjjsing < Vci consistent 
with the known presence of exponential spin correlations at finite temperatures. 

In the previous discussion we have implicitly assumed the irrelevance of terms in the free 
energy density that favor a particular surface orientation. These contributions should be 
periodic in the height and thus have a general form 

Fr=J2Vn COS 27mh (9) 
n 

where the most "dangerous" term is the case n = 1 with dimension [L]^~^. Thus for a 
stiffness constant larger than 

Vc. = \ (10) 

the surface no longer fluctuates on long length-scales but remains essentally fiat.0 

The possibility of a fluctuating defect-free surface at finite temperatures then exists for 
values of the stiffness coefficient in the region 

Vc,>V> VC2 (11) 

where defects are bound but the surface is not spontaneously oriented. In this window (|TTp 
the correlations of the discrete spins are power-law in contrast to their exponential and 
constant behavior for rj < rjc^ and rj > rjc^ respectively (see Figure 6). We note that this 
possibility for spin liquid behavior exists because the minimum Burger's vector associated 
with the underlying spin model is bmin = 6, which is twice that allowed for the pure surface 
problem. The value of the stiffness constant associated with the nearest-neighbor Ising 
model is less than the critical value required for defect-binding, consistent with the absence 
of power-law Ising spin correlations at any finite temperature. We would like to construct a 
generalization of this model with a spin-spin interaction that leads to an enhanced stiffness 
coefficient. In this paper, we have introduced additional continuous degrees of freedom and 
show that their interaction with the Ising variables leads to the stability of a classical spin 
liquid phase at finite temperatures. 



III. THE GENERALIZED SPIN MODEL 



The simplest frustrated model that has both Ising and continuous degrees of freedom is 
described by the generalized xy Hamiltonian 

^ = ^E/(0^-0.) (12) 

on a triangular lattice. We specify the periodic function /(6j — Qj) so that in the ground- 
state configurations the spins are nearly coUinear with angular orientations described by 

Gi = TTcr, + Oi (13) 

where ai and 0j are the Ising and the continuous variables respectively. The angular devia- 
tion between spin orientations is small on short length-scales (i.e. \9i — 9j \ << vr where i and 
j are nearest-neighbor sites) but can be substantial at large distances. We tune the spectrum 
of the continuous degrees of freedom to be very soft so as to maximize their interaction with 
the Ising variables. Specifically we choose the spin-spin interaction 

/(0) = ;2|cos0 + ^cos'*0|-(l-z)|cos2e-icos^0| (14) 

which admits coUinear ground-state configurations for z < z^, the coefficients have been 
selected so that the first non-vanishing term in a small-angle expansion 

/(0)~/?02 + 70' (15) 

of (|l^) is quartic (i.e. /3 = 0) while simultaneously maximizing the value of Zc- This 
specific form of the spin-spin interaction (0) leads to very fiat free-energy wells with small 
angular fiuctuations described by < 9'^ >~ . We thus expect the presence of these soft 
continuous degrees of freedom to modify the behavior of the discrete variables of the model. 
More specifically we hope that they will lead to a significant enhancement of the height 
stiffness and thus to the possibility of a spin liquid phase (Figure 6). We note that a simpler 
form of /(0) with /3 7^ in (|I3]) results in < 9^ >~ T; this leads to a numerically small 
increases of the surface stiffness that does not produce any interesting physical effects.^! 



At 2; = the ground-state of the Hamiltonian specified by ([T2D and (0) is a disordered 
nematic where the xy variables are ordered but the Ising are not. The finite z term provides 
an interaction between the discrete variables, and thus their behavior can be described by 
a fiuctuating surface as T ^ 0. For 2; 7^ there is an additional bond-bond interaction 
induced by the continuous degrees of freedom; specifically an expansion in the parameter 

results in the expression 

Afe = --^2V,a,a,(T,((ei - e,f{9k - eif)i (16) 
oi 

where {ij) and {kl) are sites belonging to the same bonds and I refers to the irreducible 



correlator. Conceptually, this interaction ([T6| ) is easier to understand as one between ferro- 
magnetic bonds; if we define 

{1 if (TjCr, = 1 
(17) 
otherwise 

then we can write 

aiajakai = ^fijfu - fij - fki - (18) 

which we insert in (|T6D. Since the total number of FM bonds is fixed in the ground-state 
manifold by the plaquette constraint, the terms linear in fij do not depend on the particular 
bond configuration and thus can be neglected. The bond-bond interaction can then be 
written as 

fbb = Aifm)fiij)f{kl) (19) 

iij){kl) 

where 

A^,,)^M) = ^m-9jne,-ei)')i (20) 

can be expressed bilinear form 

= "^gik - ^QikQii + ••■ (21) 



where 

g^,^l^{[e,-e,f) (22) 

1 2 

is independent of temperature as T — 0. Numerical evaluation of the constant A{ij)(ki) 
indicates that it is maximized by FM bonds with a common vertex that form an angle 
of and decays rapidly with increasing bond-bond separation. Therefore this induced 
bond-bond interaction ([T9| ) favors the flat (111) state which has the largest number of such 
bond configurations. The free energy associated with the entropy of the fluctuating surface 
f I sing ~ T] since ~ z^, at low temperatures it enhances the height stiffness and results 
in the sequence of transitions (paramagnet, liquid, magnet) discussed above (Figure 6). 



IV. THE CANTING TRANSITION 

At low temperatures the spins described by the generalized model of the previous section 
are energetically unstable in any collinear configuration, and the continuous variables 9i 
acquire a small (~ z) non-zero average value at each site. At T = a small-angle expansion 
of the energy about the collinear phase results in the expression 

m = E Uo^,f + he^ + ^(^J) (23) 

Jij = (1 + 3o-jaj) 

where 6ij = 9i — 9j and we have introduced an additional nematic field /i cos 20 in ([T^) to 
"control" the continuous degrees of freedom. For the 100 state each site has exactly two 
FM bonds, and the lowest eigenvalue of the coupling matrix Jij is zero. By contrast for the 
111 phase, where each site has three FM bonds, the matrix Jij has a minimum eigenvalue 
^min = —9. The lowest eigenvalue of Jij associated with the typical low-temperature state of 
the height variables is intermediate between these two extreme cases and is always negative. 

This canting instability can be illustrated by considering the angular deviation {6i) of 
a spin at any site i with three AFM bonds, where we maintain spins at all other sites in 
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their colhnear positions. From (p3D we find that the Gaussian energy associated with such 
an angular deviation is E = —z{9iY\ has a maximum at Oi = 0, and therefore the spin 
at site i is unstable to canting. A higher-order expansion of the energy p3| ) about the 
collinear phase results in an effective double-well potential associated with the continuous 
variables 9i. In the more general case where all the spins are free to cant, the angular 
deviations on different sites are correlated with a coupling determined by the matrix Jij in 
(p3|). At intermediate temperatures where an effective double- well potential develops, the 
spins are free to fluctuate between its two minima. The "freezing" of each spin in one of 
these potential wells will occur at a lower temperature, and the entropy associated with this 
transition is In 2 per site. 

The character of the transition is determined by the coupling matrix J^j which is the 
sum of two terms, one constant and the other determined by the actual conflguration of 
the surface; the latter can be viewed as the interaction between the canting angles and the 
height degrees of freedom. The uniform contribution leads to an effective (p^ theory which 
is in the Ising universality class. This flxed point behavior is stable against perturbative 
coupling between the continuous and the discrete variables which, by reasons of symmetry, 
must take the form 

F,nt (X {Vhfe^ (24) 



and is irrelevant at long wavelengths. More speciflcally, an average of (p^) over the height 
fluctuations shifts the transition temperature, and a power-counting argument applied to 
the renormalization of the ^^-interaction term 

J {{Vhof{S/Kf)dh ~ -i (25) 

indicates that it scales to zero at large distances. 

In the small-angle expansion of the energy in (p3D, the interaction between the continuous 
and the height variables is signiflcant, and thus we cannot exclude the possibility of a flrst- 
order transition or the presence of a new universality class. However numerical studies of this 
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model, described below, indicate the stability of the Ising fixed point. As a result, we note 
that the formation of the canted state has a feedback on the properties of the fiuctuating 
surface. Since the coupling matrix Jij has its lowest minimum eigenvalue for the fiat 111 
phase, the angular variables have their lowest energy in this configuration and thus modify 
the T — s> ground-state of the discrete variables. In the absence of an intervening first order 
transition, we thus expect a continuous increase in the stiffness constant rj which preserves 
the sequence of three states (paramagnet, liquid, magnet) of the Ising variables discussed 
above. In Figure 7 we show a schematic of the expected phase diagram; the relative positions 
of the canting and the liquid phase boundaries must be determined numerically. 



V. NUMERICAL DIAGNOSTICS 

In order to check these ideas we have performed Monte Carlo simulations on the model 
(p^) with z = 0.4; this is the maximum value for which no traces of ^ fluctuations were 



observed in the low-temperature spin conflgurations {zc = 0.56 for this model). Triangular 
arrays of 6^, 12^, 24^, 48^ and 96^ spins were sequentially cooled from a random conflgura- 
tion at high temperatures with 2.5 x 10^ spin flips/site per temperature point, and thermal 
averages were taken over periods of 2.5 x 10^ updates, providing 10^ approximately inde- 
pendent samples where the distribution of thermodynamic variables was studied. Some 
measurements (e.g. the surface stiffness coefficient) required longer runs (total number of 
flips per site Ntot = 10"^); in all cases the results were checked for stability with respect 
to increased number of flips/site. The speciflc heat results are displayed in Figure 8. At 
Tnem = 0.25 there exists a nematic transition into a phase with < e^*^ >^ 0, conflrmed by 
the accompanying jump in the spin stiffness shown in of Figure 9. At Tcant = 0.03 there is a 
sharp anomaly in the speciflc heat that increases with system size. This data excludes the 
possibility of an incipient flrst-order transition at low temperatures and thus we expect the 
sequence of transitions (paramagnet, liquid, magnet) described schematically in Figures 6 
and 7. 
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We associate the sharp structure in the specific heat data at low temperatures (Figure 
8) with the canting instabihty. A standard finite-size analysis (see Figure 10) confirms that 
this transition is in the Ising universality class. Briefiy for a finite system of size L that is 
smaller than the relevant correlation length ^ the specific heat scales as 

c„ ~ L^f{eL-^) (26) 

where a and v are the exponents associated with Cy and ^ in the limit L — » oo, and 
e = (^-^^) is the reduced temperature; here f{x) is a universal function. For the case of 
the 2D Ising ferromagnet a = and z/ = 1 so that ( ^61) becomes 

Cy = A\ogL + f{eL) (27) 

where 

hm f{eL) = -A\og{tL) (28) 

ensuring that the 2D Ising results are recovered in the limit of infinite system-size. In Figure 
10 we have displayed the universal function f{x) in ( pTj ) as a function of log(eL) for five 
system sizes; the resulting linear plot confirms the Ising values of a and v. Furthermore we 
vary the coefficient A in (|27|) to get the best fit to the numerical data; its optimal value 
{A = 0.50) is equal to the slope of the linear plot in Figure 9, consistent with (p8|). This 
coefficient A is the amplitude associated with the specific heat at the transition, and is 
consistent with the known value {Aising = 0.53) for the two-dimensional Ising ferromagnetic 
model.0 A crude integration of the specific heat peak at Tcant in Figure 8 leads to an 
estimated entropy S ~ 0.85 In 2; both these amplitude and the entropy results are consistent 
with approximately one Ising degree of freedom/spin at the transition. 

In order to confirm that the Tcant = 0.03 transition is due to the presence of the con- 
tinuous degrees of freedom, we apply a nematic field hnem = hcos26i to suppress their 
thermal fiuctuations. The specific heat anomaly remains sharp (Figure 11a) indicating that 
this field hnem is not conjugate to the order parameter; however it does affect the transition 
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temperature as shown in the resulting phase diagram in Figure lib. We can estimate the 
critical nematic field strength (hcrit) from the quadratic form of the energy; it is the min- 
imum value of h which maintains the expression ( |23D positive definite. For a typical state 
of Ising variables at T = 0.03, we have diagonalized the coupling matrix Jij and found that 
^min = —7.5. Thus the transition is completely suppressed if hnem > 0.5 in rough agreement 
with the numerical results (Figure 11), where we were not able to follow the phase line down 
to temperatures T < 0.015 due to problems of equilibration. 

We can also estimate the canting transition temperatures using the quartic terms in the 
small-angle expansion of the energy (^3|). For this purpose we consider again the canting 
instability of a spin at site i with three AFM bonds, where again we maintain all its neighbors 
in their coUinear positions. The renormalized energies associated with AFM and FM bonds 



are 



E^FM ^ 01 {\z + {el) 



(29) 



so that the spin at site i becomes unstable at 



where Ps{T) is the spin stiffness; ( pO]) is the defining equation for the canting transition 
temperature. Our measured quantities 

^ ^ 0-06 ^ I (31) 

are consistent with the expression (pO|). 

In order to characterize the fluctuations of the Ising degrees of freedom, we determine the 
height surface by measuring the relative height fluctuations {hi) = {{h?{r)) — {h{r))'^ between 
the middle and the edge of each array. We observe that for T < T^em the height fluctuations 
are described by the T = flxed point of the pure Ising antiferromagnet; the logarithmic 
behavior of the fluctuations, shown in Figure 12, is only robust up to the average distance 
between height defects ((^ ~ 3 x 10^) which exceeds all system sizes studied. As expected, 
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the height stiffness {^{h"^) = increases in the approach to the canting transition; this 
behavior is also displayed in the profile of the height fluctuations as a function of decreasing 
temperature (Figure 13). Unfortunately we cannot maintain the discrete degrees of freedom 
in equilibrium to low enough temperatures to directly observe the presence of the height 
liquid. However, since the specific heat data indicates no first-order transition, the surface 
stiffness must increase continuously; furthermore we know, from our arguments about the 
canting instability, that the surface will become spontaneously oriented with r] > r]c2 in 
the limit T — 0. Our numerical results confirm the expected enhancement of the surface 
stiffness coefficient with decreasing temperature. We therefore infer the presence of a finite- 
temperature fluctuating surface with no defects; its magnetic analogue is a spin liquid with 
power-law Ising correlations. 

Thus we observe that in the vicinity of the canting transition, the surface stiffness is 
enhanced. We expect that for T < Tcant V "^iH increase above its critical value (r^ci), resulting 
in a finite-temperature height liquid on all length scales. However on our length-scales of 
measurement even the Ising system is a liquid, since at T = 0.03 its correlation length is 
several thousand lattice spacing, well beyond the largest system sizes studied. We have used 
the surface stiffness to characterize the low-temperature behavior of the generalized model 
and to distinguish it from its discrete counterpart. 

The dynamics of the resulting liquid, even at short length-scales, is also very different 
from that of the pure Ising case due to the interaction between its discrete and continuous 
degrees of freedom. In order to illustrate this feature, we consider small changes to our 
fluctuating surfaces associated with the addition/subtraction of an elementary block. For 
the pure Ising case there is no energy cost associated with this process; however it changes 
the structure of the coupling matrix Jij in (]23|), and results in a dramatic reconstruction of 
the continuous degrees of freedom. As T ^ T^ant we expect that the canting instability will 
appear first in the lowest eigenmodes Jij, so that we can neglect all the others. In Figure 14 
we display a typical surface and its lowest eigenmode is shown in Figure 15. The subtraction 
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of a cube not only rearranges the energy level ordering but also modifies each state. In order 
to illustrate this property, we subtract a single block from the surface displayed in Figure 
14 at the location denoted by the arrow; we then rediagonalize the coupling matrix and 
display the eigenmode that has the maximum overlap with that shown in Figure 15 . The 
dramatic reconstruction of the continuous degrees of freedom displayed by comparison of 
Figures 15 and 16 requires some time; in the transient regime the energy of the angular 
degrees of freedom is high implying the presence of macroscopic barriers associated with 
this local process. 

The preceeding discussion suggests a distribution of surface relaxation times. @'§ We 
thus measure the time-dependent overlap 

1 ^ 

<5W = T?E(-1)'''^*^"'''^°^ (32) 
^ i=i 

between initial and final Ising states at a given temperature, probing the system's exploration 
of phase space ; any feature in Q{t) indicates memory on the time-scale of those structure^ 
Before presenting our results, let us pause briefly to discuss the possible scenarios. In a high- 
temperature paramagnet ( "gas" ) the exponential spatial and temporal spin correlations lead 
to an exponential decay of the overlap. By contrast, in a spin liquid the decay of Q{t) is 
power-law; this is the situation for the T —>■ Ising triangular antiferromagnet, and is 
maintained for this discrete case at any finite temperature up to the length-scale defined by 
the defect-defect separation. In a glass the functional form of the overlap decay is yet slower 
than power-law due to the broad distribution of relaxation times. A slow decay of the overlap 
on a time-scale r indicates the presence of barriers of height ~ Tlnr, and the magnitude of 
Q{t) measures the number of metastable states separated by these barriers. For example, 
the Sherrington-Kirkpatrick (SK) spin glass, where the interaction is infinite-range, has a 
minimum correlation volume that scales with the system-size (L); the associated overlap 
decay indicates that both the barrier number and their maximum height scale as powers of 
L. By contrast, the Edwards- Anderson spin glass with its short-ranged interaction, displays 
a broad distribution of length-scales; the associated barrier distribution is continuous and 
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only limited by L, where the number of barriers with maximum height scales as where 
a> 0. 

In Figures 17 and 18 we show our results for the time-dependent overlap of the generalized 
model above and below the canting transition; the measurements are averaged over different 
time-scales. In Figure 17 we observe that the short-time behavior of is quahtatively 

similar for T > Tcant and T < Tcant', both display similar features at these time-scales 
but are different from the pure Ising case where this is completely flat (dotted hne 

in Figure 17). At longer times (Figure 18) only the overlap measured at T < Tcant has 
discernible structure; [(3(^)]^ for T > Tcant on these time-scales is featureless. Figure 18a 
indicates that for T < Tcant the system is making transitions between a few metastable states 
separated by large barriers of similar height; their size is proportional to the logarithm of the 
pseudo-period of the observed undulations. This data also suggests that there is a minimum 
correlation volume associated with this generalized model; presumbably it is related to the 
number of sites involved in the reconstruction of the 9i variables following an elementary 
process. A study of the identical overlap for a larger system-size is necessary to determine 
the finite-size scaling of the maximum barrier height and number; in particular we would 
hke to investigate whether the system becomes "locahzed" in phase space thereby breaking 
ergodicity. The results displayed in Figure 18b indicate that the barrier heights do increase 
for larger L; however 

lim[Q(t)]2= lim[Q,,,„,(t)]2 (33) 

t—*oo t— >oo 

so that the number of metastable states does not increase with system-size and the minimum 
correlation volume remains finite. Thus, though the system retains some long-time memory 
of its initial Ising state, it is a "slow" viscous hquid and not a glass; this is because in the 
thermodynamic limit (L oo) it can explore almost all of phase space and thus is effectively 
ergodic. 
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VI. DISCUSSION 



In summary we have constructed a frustrated spin model whose discrete degrees of free- 
dom acquire power-law correlations at finite temperatures. There is a true thermodynamic 
transition into this "spin liquid", which exists as an intermediate phase separating the 
paramagnetic ("gas") and the magnetic ("sohd") states. This spin hquid has an extensive 
zero-point entropy and is the finite-temperature analogue of the T — > phase of the Ising 
antiferromagnet on the triangular lattice. This highly degenerate ground-state manifold is 
conveniently parametrized by an undulating surface, whose fluctuations on large length- 
scales are characterized by a surface stiffness constant. For the purely discrete model, the 
mapping between the spin and the membrane problem is no longer uniquely defined at any 
finite temperature due to the presence of singularities in the surface. 

In this paper we have constructed an xy spin model on the triangular lattice that has 
both continuous and Ising low-energy excitations; the ground-state configurations of these 
discrete variables are isomorphic to those of the pure Ising antiferromagnet above. In the 
generalized model, the interaction between the xy and the Ising degrees of freedom results 
in the continuous enhancement of the surface stiffness with decreasing temperature. As 
a result, the defects are bound at finite temperatures and the surface mapping becomes 
well-defined. For very low temperatures the stiffness coefficient becomes sufficiently large to 
spontaneously orient the surface, and the ground-state is unique. 

We thus expect a sequence of three phases (paramagnet, liquid, magnet) associated with 
the discrete variables as a function of decreasing temperature. There is also a canting insta- 
bility where the angular deviations 9i acquire a non-zero average at each site; the locations 
of these canting and liquid phase boundaries must be determined numerically. Similarly the 
possibility of an intervening first-order transition cannot be excluded analytically. Therefore 
we have performed Monte Carlo simulations on the generalized xy model described here. At 
low temperatures {Tcant — 0.03 for z — 0.4) we observe a sharp anomaly in the specific heat 
which we attribute to the canting transition; a finite-size scaling analysis of these results 
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indicates that it is in the Ising universahty class. Apphcation of a nematic field suppresses 
Tcant, confirming that it is due to the continuous degrees of freedom. The resulting phase 
line and the measured values of Tcant and Ps{T = Tcant) are all consistent with that expected 
from a small-angle expansion of the energy. Measurement of the surface stiffness coefficient 
confirms its enhancement with decreasing temperature, though unfortunately equilibration 
problems prevent the direct observation of the height liquid when rj > rjci. The interaction 
between the discrete and the continuous variables also leads to the presence of macroscopic 
barriers associated with local processes. In order to probe the systems exploration of phase 
space, we measure the overlap between the final and initial Ising states averaged over dif- 
ferent time-scales; we find that though the system has some long-term memory it remains 
effectively ergodic in the thermodynamic limit. 

More generally, we would like to emphasize that the classical liquid discussed here is 
not an anomalous result restricted to the specific Hamiltonian constructed above. In the 
broadest terms, the stability of such a finite-temperature phase with power-law discrete spin 
correlations could emerge from any model with the following crucial features: 

(i) At low temperatures it should have both discrete and continuous spin excita- 
tions 

(ii) Its discrete degrees of freedom retain an extensive entropyEl at T ^ and, 
in all known examples in finite-dimensions, the associated ground-state manifold 
is conveniently described by a freely fiuctuating surface 

(in) The interaction between the discrete and the continuous degrees of freedom 
should enhance the stiffness constant of the associated surface 

We note, for example, that an xy model on the kagome lattice satisfies all three of these 
conditions; furthermore since rjpous = Va we expect that a weak interaction between the 
continuous and discrete degrees of freedom should be sufficient to drive the system into a 
liquid phase. 

The enhanced stiffness constant associated with the model specified by criteria (z), {ii) 
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and (iii) not only affects the equilibrium properties of the underlying discrete degrees of 
freedom, but modifies their dynamics as well. In particular a local process, which in the 
purely discrete model has a negligible energy cost, can now involve significant rearrangement 
of the continuous degrees of freedom; this surface reconstruction is not instantaneous and 
thus there will be a transient regime where the surface energy is large. The distribution 
of the associated barriers can be probed by a study of the overlap decay between final and 
initial discrete states. For the specific model studied in this paper, we find that there is a 
large minimum correlation volume that remains finite in the thermodynamic limit resulting 
in a viscous spin liquid that effectively retains its ergodicity. Though the initial Hamiltonian 
has purely local couplings, the interplay between the Ising and the xy degrees of freedom 
generates long-range interactions between discrete spins on a length-scale that is large but 
finite in the limit L — > oo. The development of a glass requires the generation of an infinite- 
range correlation; this may be difficult to achieve in two dimensions due to the large phase 
space available to long-wavelengh spin fiuctuations. It is difficult to construct an analogous 
model in higher dimensions since the tendency towards magnetic ordering is larger than in 
d — 2, and thus it will be difficult to find a discrete system that retains its extensive entropy 
in the limit T — > and thus to satisfy (ii). However, if found, we expect it to have still 
slower dynamics and even possible glassiness in the absence of disorder. 
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FIGURES 

FIG. 1. A section of the (111) ground-state configuration of the Ising triangular model a) before 

and b) after the vr-flip of a spin (bold) with exactly three AFM bonds (solid lines); all other circled 
spins can be similarly reversed at zero energy cost resulting in a lower bound for the T — > entropy 
ofS> f ln2. 

FIG. 2. A section of the (100) ground-state configuration of the Ising triangular model a) before 
and b) after the vr-flip of a spin (bold) with four AFM bonds (solid lines); this process creates two 
"bad" triangles (shaded) whose subsequent separation by an additional spin-flip, illustrated in c), 
does not further violate the plaquette constraint. 

FIG. 3. A schematic illustrating that the energy cost associated with two defects in the Ising 
triangular model is independent of their relative position; here only the AFM bonds are displayed, 
and the dotted lines refer to the bonds "flipped" in moving a defect from its location in a) to that 
in b). 

FIG. 4. A section of a typical ground-state Ising configuration (a) with its corresponding rhombi 
tiling (b) and its associated surface representation (c) where the heights {hi) are indicated. 

FIG. 5. The height differences around a closed loop (bold line with arrows) in the a) the absence 
and b) the presence of a height defect (dotted lines indicate the bonds "flipped" to create the defect 
pair) ; note that this loop sum is zero if the number of defects enclosed is even but is a multiple of 
six if it is odd. 

FIG. 6. Schematic of the sequence of paramagnet-liquid-magnet transitions expected as a func- 
tion of increasing height stiffness constant. 

FIG. 7. Proposed phase diagram for the generalized spin model described in the text where 
C, L and S refer to the canted, liquid and solid phases repectively and the dotted line indicates a 
first-order transition . 
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FIG. 8. Specific heat results for the z = 0.4 model described in the text for arrays of 6^,12^, 
24^, 48^ and 96^ spins. 

FIG. 9. Spin Stiffness Results for the z = 0.4 model described in the text for arrays of 6^,12^, 
24^, 48^ and 96^ spins. 

FIG. 10. Finite-size Scaling Analysis for the Specific Heat Data T < Tc where = 0.03; the 
linear behavior of this plot confirms that the associated canting transition is in the Ising universality 
class (similar results are obtained for the T > Tc case). 

FIG. 11. a) Specific heat data for hnem = 0.0,0.3 indicating that the anomaly remains sharp 

despite the suppression of Tc and b) Phase diagram of hnem vs Tising where hnem is an applied 
nematic field (see text). 

FIG. 12. < h^ > vs. InL where the height fluctuations (/i^) = {h'^{r)) — {hir))"^ are related to 
the height stiffness (77) and the array length (L) by the expression < h'^ >= 

FIG. 13. Profile of height fluctuations across the diagonal of a 96 x 96 array as a function of 
decreasing temperature. 

FIG. 14. A snapshot of a "typical" surface in the vicinity of the canting transition; the arrow 
points to one location where the addition/subtraction of an elementary cube can occur. 

FIG. 15. The lowest eigenmode (|^|) of the "typical" state in Fig. 14 in the vicinity of the 
canting transition 

FIG. 16. An eigenmode corresponding to a state after the subtraction of an elementary 
building block in Fig. 14 at the location denoted by the arrow; the eigenmode displayed here has 
a maximal overlap with the original one shown in Fig. 15. 
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FIG. 17. The time-dependent height overlap [Q{t)]^ = (7^ ^.(-l){'»<(*)-'»i{o)}^ averaged over 
a) 10 and b) 100 Monte Carlo steps for a 24^ array; its behavior for T > and T < is 
qualitatively similar. 

FIG. 18. The time-dependent height overlap [Q{t)]'^ = (ly- Ei(-1)^'**^*^~''*^°^^)^ averaged over 
1000 Monte Carlo steps for ) and a b) 48^ array; the pseudo-periodic structure implies the 
presence of large barriers whose height but not number increase with system size and thus the 
system is effectively ergodic in the thermodynamic limit. 
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